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Abstract. The process of phase separation of binary systems is described 
by the Cahn-Hilliard equation. The main objective of this article is to give 
a classification on the dynamic phase transitions for binary systems using ei- 
ther the classical Cahn-Hilliard equation or the Cahn-Hilliard equation coupled 
£^ with entropy, leading to some interesting physical predictions. The analysis is 

^ based on dynamic transition theory for nonlinear systems and new classifica- 

5 tion scheme for dynamic transitions, developed recently by the authors. 
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1. Introduction 

Cahn-Hilliard equation describes the process of phase separation, by which the 
two components of a binary fluid spontaneously separate and form domains pure 
in each component. The main objective of this article is to provide a theoretical 
approach to dynamic phase transitions for binary systems. 
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Classically, phase transitions are classified by the Ehrcnfcst classification scheme, 
based on the lowest derivative of the free energy that is discontinuous at the tran- 
sition. In general, it is a difficult task to classify phase transitions of higher order, 
which appears in many equilibrium phase transition systems, such as the PVT sys- 
tem, the ferromagnetic system, superfluids as well as the binary systems studied in 
this article. 

For this purpose, a new dynamic transition theory is developed recently by the 
authors. This new theory provides an efficient tool to analyze phase transitions 
of higher order. With this theory in our disposal, a new dynamic classification 
scheme is obtained, and classifies phase transitions into three categories: Type- 
I, Type-II and Type-Ill, corresponding mathematically to continuous, jump, and 
mixed transitions, respectively; see the Appendix as well as two recent books by 
the authors [3J 0] for details. 

There have been extensive studies in the past on the dynamics of the Cahn- 
Hilliard equations. However, very little is known about the higher order transitions 
encountered for this problem, and this article gives a complete classification of the 
dynamics transitions for binary systems. The results obtained lead in particular 
to various physical predictions. First, the order of phase transitions is precisely 
determined by the sign of a nondimcnsional parameter K such that if K > 0, the 
transition is first-order with latent heat and if K < 0, the transition is second- 
order. Second, a theoretical transition diagram is derived, leading in particular 
to a prediction that there is only second-order transition for molar fraction near 
1/2. This is different from the prediction made by the classical transition diagram. 
Third, a critical length scale is derived such that no phase separation occurs at any 
temperature if the length of the container is smaller than the critical length scale. 
These physical predictions will be addressed in another article. 

This article is organized as follows. In Section 2, both the classical Cahn-Hilliard 
equation and the Cahn-Hilliard equation coupled with entropy are introduced in a 
unified fashion using a general principle for equilibrium phase transitions outlined in 
Appendix B. Sections 3-6 analyze dynamic transitions for the Cahn-Hilliard equa- 
tion in general domain, rectangular domain, with periodic boundary conditions, 
and for the Cahn-Hilliard equation coupled with entropy. Physical conclusions are 
given in Section 7, and the dynamic transition theory is recalled in Appendix A. 

2. Dynamic Phase Transition Models for Binary Systems 

Materials compounded by two components A and B, such as binary alloys, bi- 
nary solutions and polymers, are called binary systems. Sufficient cooling of a 
binary system may lead to phase separations, i.e., at the critical temperature, the 
concentrations of both components A and B with homogeneous distribution un- 
dergo changes, leading to heterogeneous distributions in space. Phase separation 
of binary systems observed will be in one of two main ways. The first is by nuclc- 
ation in which sufficiently large nuclei of the second phase appear randomly and 
grow, and this corresponds to Type-II phase transitions. The second is by spinodal 
decomposition in which the systems appear to nuclear at once, and periodic or 
semi-periodic structure is seen, and this corresponds to Type-I phase transitions. 

Since binary systems are conserved, the equations describing the Hclmholtz pro- 
cess and the Gibbs process are the same. Hence, without distinction we use the 
term "free energy" to discuss this problem. 



CAHN-HILLIARD EQUATIONS AND PHASE TRANSITION DYNAMICS FOR BINARY SYSTEMS 



Let ua and Ms be the concentrations of components A and B respectively, then 
ug = 1 — ua. In a homogeneous state, ug = ug is a constant, and the entropy 
density So = Sq is also a constant. We take u, S the concentration and entropy 
density deviations: 

u = ub~u Bi S=S -S 
By ( B. 1 1 and (B.2|, the free energy is given by 
(2.1) 

F{u, S) = F Q - 



^ 'V^l 2 + f |VSf + f |Sf + (3 Su + (3 2 Su 2 + a k u k 



2p 



fc=l 



dx. 



Since entropy is increasing as u —* 0, and by 

-F(u, S) = -^ 2 AS + p x S + o u + (3 2 u 2 = 0, 



5S 



we have 
(2.2) 



A, = 0, /3i > 0, (3 2 > 0, 



which implies that S is a decreasing function of \u\ 



According (B.lll and (B.5l, we derive from (2.1) and (2.2) the following equa- 



tions governing a binary system 

as 
at 

du 



(2.3) 



dt 



= kiAS — aiS — a 2 u 2 , 
= -k 2 A 2 u + b A(Su) + Af(u), 



u{x, t)dx = 0, 



where fci, k 2 , bo, a\, a 2 are positive constants, and 

2p-l 

(2.4) f(u) = ^ M fc , 6 2p _! > (p > 2). 



fe=i 



Physically sound boundary condition for (2.3 1 is either the Neumann boundary 
condition: 



= 0, — = on oil, 



dn dn ' dn 
with il C K™ (1 < n < 3) being a bounded domain, or the periodic boundary 
condition: 

(2.6) u(x + KL) =u{x), S(x + KL) = S{x) 

with ft = [0,L] n ,K= (fci,--- ,fc„),l <n<3. 

For simplicity, in this section we always assume that p — 2. Thus, function (2.4) 
is rewritten as 

(2.7) f(u) = b lU + b 2 u 2 + b 3 u 3 , b 3 >0. 

Based on Theorem | A. 1[ we have to assume that there exists a temperature T\ > 
such that b\ = b\(T) satisfies 

>0 HT>T U 

(2.8) b x (T) \ =0 if T = T X , 

<0 if T<T X . 
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If we ignore the coupled action of entropy in (2.1 1, then the free energy F is in 
the following form 



F(u) =F + 



Ml 



Vul 



"1 2 , "2 3 Q!3 4 



da;, 



and the equations (2.3) are the following classical Cahn-Hilliard equation: 

du 
dt 



(2.9) 



■kA 2 u + Af(u), 
u(x, t)dx = 0, 



where f(u) is as in (2.7 1. 

3. Phase Transition in General Domains 

In this section, we shall discuss the Cahn-Hilliard equation from the mathemat- 
ical point of view. We start with the nondimensional form of equation. Let 



I 4 , 



X 



l 2 h 



72 



l 2 b 2 u 



U = UqU , 

73 ~ 



k k k 

where I is a given length, uq = ug is the constant concentration of B, and 73 > 0. 
Then the equation (2.9 1 can be rewritten as follows (omitting the primes) 

-A 2 u - XAu + A(j 2 u 2 + 73U 3 ), 
u(x, t)dx = 0, 



du 
dt 



(3.1) 



Let 



u(x, 0) = ip. 

H = |« S L 2 (fl) I J udx 
iary condition ( 

Hi = ju G i? 4 (Sl) n H 



For the Neumann boundary condition (2.5) we define 

du dAu 



dn dn 



on dfl 



and for the periodic boundary condition (2.6) we define 

= H 4 ((l) C\H u(x + KL) = u{x) MK G Z n | . 

Then we define the operators L\ = —A + B\ and G : H\ — > i? by 

Aw = A 2 u, 

(3.2) S A ?i = -AAtt, 

G(u) = 72 Au 2 +73AU 3 . 



Thus, the Cahn-Hilliard equation (3.1) is equivalent to the following operator 
equation 

du T 

L\U + Cr(U), 



(3.3) 



dt 

u(0) = <fi 
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It is known that the operators defined by (|3.2| satisfy the conditions (A. 2) and 



(A.3I 



We first consider the case where C E™(1 < n < 3) is a general bounded 
and smooth domain. Let pk and be the eigenvalues and eigenfunctions of the 
following eigenvalue problem: 

- Ae k = Pkdk, 
de k , 



(3.4) 



dn 



I an 



0, 



efedx = 0. 



The eigenvalues of (|3.4| satisfy < pi < p 2 < • • ■ < Pk < • • • , and lim fe ^ 

oo Pk — OO- 



The eigenfunctions {e^} of (3.4 1 constitute an orthonormal basis of H. Further- 
more, the eigenfunctions of (3.4 1 satisfy 



dAe k , 
dn 



\dn 



0, fc = l,2, 



Hence, {e^} is also an orthogonal basis of Hi under the following equivalent norm 

-l 1/2 



Mi = 



lA^urdx 



We are now in position to give a phase transition theorem for the problem (3.1 ) 
with the following Neumann boundary condition: 



(3.5) 



du 
On 



= 0, 



dAv 
dn 



= 



on dCl. 



Theorem 3.1. Assume that 72 
tions hold true: 



and 73 > in (3.1), then the following asser- 



(1) If the first eigenvalue p\ of (3.4-) has multiplicity m > 1, then the problem 
(3.1) with (3.5) bifurcates from (u, A) = (0,pi) on A > p\ to an attractor 
Yi\, homeomorphic to an (m — Y)- dimensional sphere S m , andT,\ attracts 
H \ T, where T is the stable manifold of u = with codimension m. 

(2) T,\ contains at least 2m singular points. If m — 1, has exactly two steady 
states ±U\, and if m = 2, Y,\ = S 1 has at most eight singular points. 

(3) Each singular point u\ in £>, can be expressed as 

u x = (X~p 1 ) 1/2 w + o(\X- Pl \ 1 / 2 ), 



where w is an eigenfunction corresponding to the first eigenvalue of (3.4) 
Proof. We proceed in several steps as follows. 



Step 1. It is clear that the eigenfunction {e^} of (3.4) are also eigenvectors of 
the linear operator L\ = —A + B\ defined by (3.2 1 and the eigenvalues of L\ are 
given by 



(3.6) 



Ac(A) = p fe (A~p fc ), fc=l,2,---. 



It is easy to verify the conditions (A.4| and (A. 5 1 in our case at A = p\. We shall 
prove this theorem using the attractor bifurcation theory introduced in [3]. 
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We need to verify that u — is a global asymptotically stable singular point of 



(3.3| at A = pi- By 72 = 0, from the energy integration of (3.1 ) we can obtain 

u 2 dx 



(3.7) 



1 d 
2dt 



[ u 2 dx = [ [-|Aw| 2 + p 1 |Vw| 2 - 37 3 u 2 |Vw| 2 ] dx 
Jn Jn 



< -C 



[ |Av| 2 da;- 373 / u 2 \\7u\ 2 dx, 
Jn Jn 



where C > is a constant, u = v + w, and vwdx — 0, and w is a first cngcnfunc- 
tion. It follows from (3.7) that u — is global asymptotically stable. Hence, for 
Assertion (1), we only have to prove that Sa is homeomorphic to S™ -1 , as the rest 
of this assertion follows directly from the attractor bifurcation theory introduced 
in 0. 



Step 2. Now we prove that the bifurcated attractor from (0, pi) contains at 
least 2m singular points. 

Let g(u) = —Ait — Xu + 73U 3 . Then the stationary equation of (3.1 1 is given by 



Ag(u) = 0, 
udx = 0, 



which is equivalent, by the maximum principle, to 



(3. 



Au — Xu + 73 it = constant, 



udx = 0, 



m 
du , 
On 



I dn 



0. 



By the Lagrange multiplier theorem, (3.8 1 is the Euler equation of the following 
functional with zero average constraint: 



(3.9) 



1 



iVul- 



F(u)-- 



X 



73 4 



dx. 



du , 
dn 



\on 







udx = 



Since F is an even functional, by the classical Krasnoselskii bifurcation theorem 
for even functionals, (|3.9| bifurcates from A > p\ at least to 2m mini-maximum 



points, i.e., equation (3.8 1 has at least 2m bifurcated solutions on A > p\. Hence, 



the attractor contains at least 2m singular points. 



Step 3. To complete the proof, we reduce the equation (3.3 1 to the center 
manifold near X — p\. By the approximation formula given info], the reduced 
equation of (3.3 1 is given by: 



(3.10) 



dx 



l3 1 {X)x l - 73P1 / v 3 eudx + o(|x| 3 ) 
at J n 



for 1 < i < m. 
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where ft (A) = pi(A - pi), v = YT =-\ an d {en, • • • 

functions of (3.4). Equations (3.10[) can be rewritten as 



, ei, 



t } are the first eigen- 



(3.11) 



Let 



dxi 
~dt 



= (3l(X)Xi - 73/91 




eudx + o(|a;| 3 ). 



g(x) = I J v eudx, 
Then for any x £ R m , 
(3.12) <g(x),x> 



v 3 ei r 



, dx 



v(x) 



3=1 



Cjeij. 



8 = 1 



v 3 eudx — / w 4 cfa: > C|x| 4 , 



for some constant C > 0. 

Thus by the attractor bifurcation theorem [3], it follows from (3.111 and (3.12) 
that the attractor T,\ is homeomorphic to S™ -1 . Hence, Assertion (1) is proved. 
The other conclusions in Assertions (2) and (3) can be derived from (3.111 and 
(3.121. The proof is complete. □ □ 

Physically, the coefficients 72 and 73 depend on uq — ub, the temperature T, 
and the pressure p: 

7*, =7fc(uo,T,p), fc=2,3. 
The set of points satisfying 72(^0, T,p) = has measure zero in (uo,T,p) £ R 3 . 
Hence, it is more interesting to consider the case where 72 7^ 0. 

For this purpose, let the multiplicity of the first eigenvalue pi of (3.4 1 be m > 1, 
and {ei, • • • , e m } be the first eigenfunctions. We introduce the following quadratic 
equations 



(3.13) 



Theorem 3.2. 



ft; 







for 1 < k < m, 



etejekdx. 
Let 72 7^ 0, 73 > 0, and x 



(3.13). Then the phase transition of (3.1 
Furthermore, the problem (3.1) with 



be an isolated singular point of 
) and (3.5) is either Type-II or Type-Ill. 
3.5) bifurcates to at least one singular point 



on each side of X — pi, and has a saddle-node bifurcation on A < p\. In particular, 
if m = 1, then the following assertions hold true: 

(1) The phase transition is Type-Ill, and a neighborhood U C H of u — can 
be decomposed into two sectorial regions U — Di(n) + such that the 
phase transition in Di(ir) is the first order, and in .02(71") is the n-th order 
with n > 3. 

(2) The bifurcated singular point u\ on A > p\ attracts Z3 2 ( 7r ), which can be 
expressed as 

(3.14) u\ = (A — pi)ei/7 2 ft + o(|A - Pl \), 

where, by assumption, a = J„ efdx 7^ 0. 

(3) When I720I = s is small, the assertions in the transition perturbation theo- 
rems (Theorems A. 8 and\A.(fy hold true. 
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Remark 3.1. We shall sec later that when fl is a rectangular domain, i.e., 

ft = n£ =1 (0, Lk) C K™ (l<n<3), 

then a = Jq efdx = 0. However, for almost all non-rectangular domains f2, the first 
eigenvalues are simple and 0. Hence, the Type-Ill phase transitions for general 
domains are generic. □ 



Proof of Theorem 3.2 Assertions (l)-(3) can be directly proved using Theorems 
A. 5 |A.8| and A. 9 By assumption, u = is a second order non-degenerate singular 



point of (3.3 I at A — pi, which implies that u — is not locally asymptotically sta- 
ble. Hence, it follows from Theorems |A.3| and the steady state bifurcation theorem 
for even-order nondegenerate singular points [3] that the phase transition of (3.1) 
with (3.5 1 is either Type-II or Type-Ill, and there is at least one singular point 
bifurcated on each side of A = p\ . 

Finally, we shall apply Theorem |A.6| to prove that there exists a saddle- node 
bifurcation on A < p\ . 

It is known that 

{even if A = px, 
1 ifA<p, 



Moreover, since L\ + G defined by (3.2 1 is a gradient-type operator, we can derive 
that 

md(L pl + G, 0) < 0. 
Hence there is a bifurcated branch £ A on A < p\ such that 

ind(L A + G, u\) = -1 Vw A e S A , A < p x . 



It is clear that the eigenvalues (3.6) of L\ satisfy (A.4| and (A.5I. For any A £ K 
(3.3 1 possesses a global attractor. Therefore, for bounded A, 6 < A < pi, the 
bifurcated branch S A is bounded: 



\u\\\h<C Vua £ £; 



-oo < b < A < px. 



We need to prove that there exists A < pi such that for all A < A equation (3.3) 
has no nonzero singular points. 

-72/273 and u ^ in H, 



By the energy estimates of (3.3), for any A < A 
[|Au| 2 -A|Vu 

\Au\ 2 dx 



j 

Jo. 



27 2 u|Vw| 2 + 37 3 u 2 |Vu| 2 ] dx 

[ |Vu| 2 (-A-2| 72 u| + 373W 2 )dx 
Jn 



> 



Wurdx- 



|V U | 2 (-A + 7 3 ^ + 2 73 ( U - ^) 

73 



± 
273 



)dx 



> 
> 



f \Vu\ 2 dx+ [ |Vu| 2 (-A 
Jn Jn 



± 
273 



)dx 



Therefore, when A < A, (3.3) has no nontrivial singular points in H. Thus we infer 
from Theorem |A.6| that there exists a saddle-node bifurcation on A < px . This 
proof is complete. □ 
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Figure 4.1. Type-II 
transition as given by 
Theorem 4.1 



4. Phase Transition in Rectangular Domains 

The dynamical properties of phase separation of a binary system in a rectangular 
container is very different from that in a general container. We see in the previous 
section that the phase transitions in general domains are Type-Ill, and we shall 
show in the following that the phase transitions in rectangular domains are cither 
Type-I or Type-II, which are distinguished by a critical size of the domains. 

Let ft = n£ =1 (0, L k ) C E™ (1 < n < 3) be a rectangular domain. We first 
consider the case where 



(4.1) 



L = Li > Lj V2 < j < n. 



Theorem 4.1. Let f2 = II£ =1 (0, L k ) satisfy (4-D- The following assertions hold 
true: 

(1) If 

2L 2 2 

73 < g^, 



then the phase transition of (3.1) and j 3.5) at A = Ao = tt 2 /L 2 is Type-LL. Ln 
particular, the problem (3.1) with (3.5) bifurcates from (u, A) = (0,7r 2 /L 2 ) 
on A < 7r 2 /L 2 to exactly two equilibrium points which are saddles, and there 



are two saddle-node bifurcations on A < tt 2 /L 2 as shown in Figure 4-1 

(2) // 

2L 2 2 

73 > 9^ 72 ' 

then the transition is Type-L. Ln particular, the problem bifurcates on A > 
7r 2 / L 2 to exactly two attractors uf and uj which can be expressed as 



(4.2) 



'1.2 



± 



V^A-f 2 ) 1 / 2 



r l/2 



TTXi 

17 



o(|A 



L 2 



whe 



_ 373 



3tt 2 



Proof. With the spatial domain as given, the first eigenvalue and eigenfunction of 

TTXi 



(3.4 1 are given by 



Pi = ttV^ 



ei = cos 



L 



The eigenvalues and eigenfunctions of L\ = — A+B\ defined by (3.2) are as follows: 

f3 K = \K\ 2 (X-\K\ 2 ), 

klTTXl k^TTX^ 



(4.3) 
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where 

K=( k f 1 --- k -f-\ Vfc.gZ, l<i<n, 

n 
i=l 

By the approximation of the center manifold obtained in [3], the reduced equation 
of (3.3 1 to the center manifold is given by 

dy 2tt 2 f 
(4.5) -£ = /3i(A)t/-— — — = / [72y 2 e? + 732/ 3 ef + 72(ye 1 +$(y)) 2 ei]cis ) 



dt (Lx---L n )L\ 
where y 6 R 1 , <&(?/) is the center manifold function, and 



(4.6) /MA) = -^(A--^). 



Direct calculation implies that 



(4.7) / e 3 ^ = • • • f L " cos 3 : : 0. 

Jn Jo Jo ' 



L 

■ L 



(4.8) / e\dx = L 2 ---L n [ cos 4 '^p-dx 1 = \l x ■ ■ ■ L n . 
Ja Jo L ° 

By pj] ) and <J>(y) = 0(\y\ 2 ) we obtain 

(4.9) / (ye 1 +<S>) 2 e 1 dx = 2y [ <S>{y)e\dx + o(\y\ 3 ) 
Jn Jn 

It follows that 

oo 

Hv)= ( t ) K(y)e K + o(\y\ 2 ) 

\K\ 2 >7T 2 /L 2 

I2V 2 f A 2 , \K\ 2 -y 2 y 2 



K ^ = r '\L 112 / Ae i • e K dx = Pi] ,Z M2 / amldx. 



Notice that 



exe 2 dx 



Li---L„/4 VX= (^,0,--- ,0 



Then we have 



Inserting ^(y) into (4.9 1 we find 

„3 



(4.10) / (yei + <S>) 2 ei dx / cos ^e?dx + o(\y\ 3 ) 

Jn A — rs- Jn L, 

y 3 + o(\y\ 3 ). 



-J2- Jn 



4 X 47T " 2 
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Finally, by (4.7 1 , (4.8 > and (4.101, we derive from (4.51 the following reduced equa- 



tion of (3.3 1 
(4.11) 



dy 
dt 




y 3 + o(\y\ 3 ). 



Near the critical point Ao = tt 2 /L 2 , the coefficient 

373 , ll _ 373 i 2 7l 



A 



4tt 2 
L 2 



3tt 2 



Thus, by Theorem A. 2 we derive from (4.11 1 the assertions of the theorem except 



the claim for the saddle-node bifurcation in Assertion (1), which can be proved in 
the same fashion as used in Theorem |3.2| The proof is complete. □ 



We now consider the case where f2 = nj? =1 (0, Lk) satisfies that 



(4.12) 



L = Li = ■ ■ ■ = L m > Lj for 2 < m < 3, m < j < n. 



Theorem 4.2. Let = n^ =1 (0, Lk) satisfy (4-. 12). Then the following assertions 
hold true: 

(1) V 



73 > 



26L 2 
27tt 2 



72 . 



then the phase transition of the problem (3.1) with (3.5) at Ao = n /L is 



Type-I, satisfying the following properties: 

(a) The problem bifurcates on A > n 2 / L 2 to an attractor containing 
exactly 3 m — 1 non-degenerate singular points, andY,\ is homeomorphic 
to an (to — 1)- dimensional sphere S" 1 ^ 1 . 

(b) For to = 2, the attractor = S 1 contains 4 minimal attractors, as 
shown in Figure \4~M 

(c) For to = 3,Sa = S 2 contains 8 minimal attractors as shown in Fig- 
ure 4-3^a), if 

22L 2 

73 < 



2 72, 



9tt 

and contains 6 minimal attractors as shown in Figure \4~3^ b ) if 

22L 2 



73 > 



(2) // 



73 < 



26L 2 
27^r 



r72> 



then the transition is Type-II. In particular, the problem has a saddle-node 
bifurcation on X < Xq = ir 2 /L 2 , and bifurcates on both side of X = Ao to 
exactly 3'™ — 1 singular points which are non-degenerate. 

Proof. We proceed in several steps as follows. 

Step 1. Consider the center manifold reduction. It is known that the eigenvalues 
and eigenfunctions of L\ = —A + B\ are given by (4.3) and (4.4 1 with L\ — ■ ■ ■ — 
L m . As before, the reduced equations of (3.3 1 are given by 

dy 



(4.13) 



dt 



^PiWy + 9{y) + o{\y\ i ), 
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where y — (yi, ■ ■ ■ ,y m ) S M. m , (3i(X) is as in (4.6 1, and 



2tt 2 



L\ ■ ■ ■ L n L 1 



(G 2 (y) + G 3 (y) + G 23 (y)), 



G 2 (y) = -72 / v exdx, ■ ■ ■ , / v e m dx , 
\Jn Jn J 

Gz{y) = -73 / v^eidx, ■ ■ ■ , / v 3 e m dx I , 
\Jn Jn J 



G 2 j,(y) = -72 ( / u eidx, ■ ■ ■ , / u e m dx 
In Jn 
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(4.14) 



v e;dx 



Here = cosirxi/L for 1 < i < m, L is given by (4.121, v — X)I=i Vi e %i u = v + ^iv) 
and <f> is the center manifold function. Direct computation shows that 

v 2 eidx= I y yj cos itXj / L I cos nXi/Ldx = 0, 

/ i/fc cos nxj I L cos 7TXi / Ldx 

Jn \j=i J 

= 2VV / $( 2 /)e i e i da;+ / <J> 2 (y) ei cfa;. 
„•_ -, Ai Jn 



3 

T 



(4.15) 



We need to compute the center manifold function <&(?/). As in [5], we have 
(4.16) 



HV)= E <My)eA' + o(|y| 2 ) + 0(| 2 /| 2 |/3 1 | 

\K\>ir 2 /L 2 



where 



(y) 



12 



-PkW < eK,e K > 
(3kW < e K ,e K > 

\K?12 



Av exdx 



v 2 exdx 



y2 y%yi I e l e i e K dx. 



It is clear that 



if 



TTXi TTXj 



eie^exdx = I cos — — cos — — e^dx = 



K^Ki + K 



1 1 



Ki 



1^1 l^n 



By (4.3) and (4.4) we have 

<f>K{y) = 



72^j(2 - Sjj) 
(\-\K\i) <e K ,e K > 

722/. 



eiejexdx 



/2j' 2 

2(A - TJ 



4^ i£ i=j, 

L 2 

2 722/i^ 

27T 2 

L 2 
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for K = Ki + Kj. Thus, by (4.16 1, we obtain 



m = y — n— y ^ cos f^i + y 



272 



„ 2v 2.yiy r cos — cos 

l>r TT> 



TTXl TTX r 



Inserting <&(?/) into (4.151 we derive 

72 



u eidx = 



cos — cos — — + 4y 4 2^ V 3 cos cos — 

3¥* 



dx 



+ o(\y\ 3 ). 
Direct computation gives that 

72L1 ■ • • L T 



(4.17) 



u eidx 



4y 4 



i 2 



+ «(M 3 )- 



Putting (4.141 and (4.17) in (4.13), we get the reduced equations in the following 
form 



(4.18) 

where 
(4.19) 



^ a (\\ 71-2 



^iy 3 + ^Vi^Vj 



373 



7l 



A 



4tt 2 



<7 2 = 373 



+ o(|y| 3 ) Vl<i<m, 



47l 

\ 2tt 2 
L 2 " 



Step 2. It is known that the transition type of (3.3 1 at the critical point Ao = 
tt 2 /L 2 is completely determined by (4.181, i.e., by the following equations 



(4.20) 
where 



dy l 
dt 



2L2 



VI < i < m, 



373 ^ 2 7 2 2 
2 3^ 2 



a° 2 = 3 73 - 



4L 2 7 2 2 



It is easy to see that 
(4.21) 



4 + 4 > 73 > , 
a? + ^ < o 73 < ^2 72- 



Step 3. We consider the case where m = 2. Thus, the transition type of (4.20) 
is equivalent to that of the following equations 

dyi 
dt 



(4.22) 



-yiKyl + ^yl] 



dy 2 
dt 



We can see that on the straight lines 
(4.23) 



2 2 
Vi = %> 
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equations (4.221 satisfy that 



dv2 
dyi 



Vi 



for irJ + o^O, (»i,|fc)^Q. 



Hence the straight lines (4.23) are orbits of (4.221 if crj + u\ 0. Obviously, the 
straight lines 



(4.24) 



?/i=0 and y 2 = 



are also orbits of (4.221 



There arc four straight lines determined by (4.231 and (4.24), and each of them 



contains two orbits. Hence, the system (4.221 has at least eight straight line orbits 



Hence it is not hard to see that the number of straight line orbits of (4.22), if finite, 
is eight . 



Since (3.3) is a gradient-type equation, there are no ellipti c regi ons at y = 0; see 
j3]. Hence, when <r° + (jj > all the straight line orbits on (4.23) and (4.241 tend 
to y — 0, as shown in Figure 4.4 (a), which implies that the regions are parabolic 
and stable, therefore y = is asymptotically stable for (422). Accordingly, by 
the attractor bifurcation theorem in [3J, the transition of (4.18) at Ao = ir 2 /L 2 is 
Type-I. 



When a I + <r\ < and a® > 0, namely 



2L 2 7 f 



9 ir< 



< 73 < 



26 

27 7T 2 ■ 



the four straight line orbits on (4.231 are outward from y = 0, and other four 
on ( 4.24| ) are toward y = 0, as shown in Figure |4.4| (b) , which implies that all 



regions at y — are hyperbolic. Hence, by Theorem |A.3| the transition of (|4.18|) 
at A = tt 2 /L 2 is Type-II 



When crj < 0, then a\ < too. In this case, no orbits of (4.22 ) are toward y = 0, 
as shown in Figure 4.4 (c), which implies by Theorem A. 3 that the transition is 
Type-II. 



and Type-II if 73 < 2 ^- 2 
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Step 4. Consider the case where m = 3. Thus, (4.20) are written as 

dm 

dt 

(4.25) 



-yi [*?!/?■ 



dt 

1 1 is clear that the straight lines 
(4.26) W = 0, = 



d V2 r 2 

= -2/2Fll/2 



^°(y 2 
^°(y 2 



% 2 )L 
% 2 )L 

2/2 2 )]- 



for i ^ j, 1 < i, j < 3, 



(4.27) 



y? = Vj, Vk = for i ^ j, i ^ k, j ^ k, 1 <i,j,k< 3, 



y 2 



?y 2 2 



y 2 . 



consist of orbits of (4.25 1. There are total 13 straight lines in (4.26 ) and (4.27 1, each 
of which consists of two orbits. Thus, (4.25) has at least 26 straight line orbits. We 
shall show that (4.251 has just the straight line orbits given by (4.261 and (4.271. 
In fact, we assume that the line 

y 2 = ZiVi, V3 = z 2Vi (zi,z 2 are real numbers) 



is a straight line orbit of (4.25). Then z±,Z2 satisfy 

! + <J 2 >(l + Z 2 ) 



Zi = zi 



(4.28) 



dm 
dyi 

dy3 t 

dm 22 Z2 ^° 1+ aUz! + zl) 



aU4{4 + zl) ' 

+ ( T»(1+Z 1 2 ) 



It is easy to see that when a® 7^ <J% the solutions z\ and Z2 of (4.281 take only the 
values 

0,±1; z 2 = 0,±l. 



z\ 



In the same fashion, we can prove that the straight line orbits of (4.251 given by 

yi = aiy 3 , y2 = c*2y3, and y x = /3ij/ 2 , y3 = fhV2 
have to satisfy that 

a 4 = 0,±l and /3 f = 0,d=l (i = 1,2). 



Thus, we prove that when <7° ^ a®, the number of straight line orbits of (4.25) is 
exactly 26. 

When a 1 = a®, we have that 73 — 7^^2-7! which implies that cr° = 0% > 0. In 
this case, it is clear that y — is an asymptotically stable singular point of (4.251. 
Hence, the transition of (4.181 at Ao = ir 2 /L 2 is I-type. 

When a^+a^ > and a { { ^ a^, all straight line orbits of (4.25 1 are toward y = 0, 
which implies that the regions at y = 0, are stable, and y = is asymptotically 
stable; see [3]- Thereby the transition of (4.181 is Type-I. 

When a® 
regions at y 



u\ < with crj > 0, we can see, as in the case of m — 2, that the 
are hyperbolic, and when a° + 0% < with u\ < the regions at 
y = are unstable. Hence, the transition is Type-II. 



Step 5. We prove Assertion (1). By Steps 3 and 4, if 73 > f^r72, the reduced 
equation (4.18) bifurcates on A > A = ir 2 /L 2 to an attractor All bifurcated 
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equilibrium points of (3.3 1 are one to one correspondence to the bifurcated singular 
points of (4.181. Therefore, we only have to consider the stationary equations: 

2 

(4.29) Pi(\)y t - ^[viyf + a 2 y l J2y';]+o(\y\ 3 ) = 0, 1 < % < m, 



where o\ and a 2 are as in (4.191. 

Consider the following approximative equations of (4.29) 

(4.30) Px(X)Vi ~ Viiaiyf + a 2 Y, Vj) = for 1 < i < m, 

where cti = 7r 2 eri/2L 2 , a 2 — ir 2 a 2 /2L 2 . It is clear that each regular bifurcated 
solution of (4.30) corresponds to a regular bifurcated solution of (4.291. 

We first prove th at (|4.30 1 has 3 m — 1 bifurcated solutions on A > Ao- For each 
k (0 < k < m - 1), ( |4.30| has C* x 2 m - k solutions as follows: 

Vn = 0, • • • , y jk = Q for 1 < ji < m, 

( 4 - 31 ) 2 2 1A A _i 

2£i = ••• = y r ^ ■ 



0i(ax + (m - k - l)a 2 ) for r,- ^ j;. 



Hence, the number of all bifurcated solutions of (4.30) is 

m— 1 

2 Ct x 2 m ~ fe = ( 2 + !)" - 1 = 3 m - 1. 



fe=0 



We need to prove that all bifurcated solutions of (4.30 1 are regular. The Jacobian 



matrix of (4.30) is given by 

fPx~hx{y) 2a 2 y 1 y 2 
2a 2 y 2 y x px - h 2 (y) 



(4.32) 



where 



Dv 



\ 2a 2 y m yi 2a 2 y m y 2 

hi(y) = 3aiy 2 + a 2 ^y 



2a 2 y x y m ^ 
2a 2 y 2 y m 



Pi - h m (y)J 



For the solutions in (4.31 1, without loss of generality, we take 

2/0 = {vi,--- ,y„), 

Vi = 1 < i < k, 

Vk+i = ■■■ = y° m = /?i 1/2 («i + (m- k - l)a 2 r 1/2 
Inserting them into (4.321 we find 

(4.33) Dv(y ) = 

where 
p = Pi ( 1 



(m — fc)a2 



Ph 
A m _ fe 

(ui - a 2 )Px 



ax + (m — k — l)a 2 J a± + (m — fc — 1)02 ' 
//3i-(3ai + (m-fc-l)a 2 )(^ +1 ) 2 ••• 



4 , — 



V 



2a 2 (^ +1 ) 2 



2a 2 (y° k+1 ) 2 
Px-(3ax + (m-k-l)a 2 )(y° k+1 ) 
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Direct computation shows that 



detA 



m— k 



(ai + (m - k - l)a 2 ) m L 2m 



det 



/-<7i cr 2 

f 2 — Cl 



\ C 2 (72 



(72 ^ 
<72 



where ai,(T2 are given by (4.19) 
re on] 

/3(A) 



Obviously, there are only finite number of A > ir 2 /L 2 satisfying 

(<7 1 (A)-<7 2 (A))/3 1 (A) 



<7 1 (A) + (m-k- 1)(7 2 (A) 
detA m _ fe (A) = 0. 



0. 



Hence, for any A — 7T /L > sufficiently small the Jacobian matrices (4.321 at the 



singular points (4.311 are non-degenerate. Thus, the bifurcated solutions of (4.30) 
are regular. 



Since all bifurcated singular points of (3.1l with (3.5 1 are non-degenerate, and 
when is restricted on XiX 3 -plane (1 < i,j < m) the singular points are con- 
nected by their stable and unstable manifolds. Hence all singular points in Y,\ are 
connected by their stable and unstable manifolds. Therefore, £a must be homeo- 
morphic to a sphere S" 1 " 1 . 

Assertion (1) is proved. 

Step 6. Proof of Assertions (2) and (3). When m = 2, by Step 5, E A = S 1 
contains 8 non-degenerate singular points. By a theorem on minimal attractors in 
[3], 4 singular points must be attractors and the others are repellors, as shown in 
Figure [O] 

When m — 3, we take the six singular points 

±Yi - (±Pia^\ 0, 0), ±Y 2 = (0, ±Pia^\ Q),±Y 3 = (0, 0, ±/3ia^) 



Then the Jacobian matrix (4.32) at Yi (1 < % < 3) is 



Pi 
Dv{±Y l ) = \ p 2 

P3 

where ; = (3\ (\ — as j ' ^ i and ^ = — 2f3\. Obviously, as a 2 < <J\, < pj 
(j ^ i) and pi < 0, in this case, ±Yfc (1 < fe < 3) are repellors in = S 2 , which 



implies that Sa contains 8 attractors ±Zk (1 < k < 4) as shown in Figure 4.3 
(a). As g 2 > <Ji,pj < (1 < j < 3), the six singular point ±Y k (1 < k < 3) are 
attractors, which implies that Sa contains only six minimal attractors as shown in 



Figure 4.3 (b). Thus Assertion (2) is proved. 

The claim for the saddle-node bifurcation in Assertion (3) can be proved by using 
the same method as in the proof of Theorem |3.2| and the claim for the singular 
point bifurcation can be proved by the same fashion as used in Step 5. 

The proof of this theorem is complete. □ 

Remark 4.1. For the domain = [0, L] m x D C R"(l < m < n) w here n > 2 is 
arbitrary and D C K n_m a bounded open set, Theorems 4.1 and 4.2 are also valid 
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provided ir 2 /L 2 < Ai, where Ai is the first eigenvalue of the equation 

- Ae = Ae, xeDc R n ""\ 
de 



dn 



OD = 0, 



edx = 0. 

D 

Remark 4.2. In Theorem |4.2[ the minimal attractors in the bifurcated attractor 
Yi\ can be expressed as 

.2 / r 2\l/2 



(4.34) u\ = (X— n'/L J ) x ^e + o(|A - ^ 2 /i 2 | 1/2 ), 

where e is a first eigenfunction of (3.4). The expression ( 4.34| can be derived from 
the reduced equations (4.181. 



We address here that the exponent (3 — 1/2 in (4.341, called the critical exponent 
in physics, is an important index in the phase transition theory in statistical physics, 
which arises only in the Type-I or the continuous phase transitions. It is interesting 



to point out that the critical exponent /3 = 1 in (3.141 is different from these (i = 1/2 



appearing in (4.2 1 and (4.341. The first one occurs when the container C K is a 



non rectangular region, and the second one occurs when f2 is a rectangle or a cube. 
We shall continue to discuss this problem later from the physical viewpoint. □ 

5. Phase Transitions Under Periodic Boundary Conditions 

When the sample or container £1 is a loop, or a torus, or bulk in size, then the 
periodic boundary conditions are necessary. In this section, we shall discuss the 
problems in a loop domain and in the whole space = W 1 . 

Let ft — S 1 x (n, r 2 ) C i 2 be a loop domain, < r% < r 2 . Then the boundary 
condition is given by 

u(9 + 2/ctt, r) = u(9, r) for < 6 < 2ir, r 1 <r < r 2 , fee Z, 

(5-1) du d 3 u 

at r — ri,r 2 - 



Or 



or 6 



Assume that the gap r 2 — r\ is small in comparison with the mean radius tq = 
(Vi + r 2 )/2. With proper scaling, we take the gap r 2 — r\ as r 2 — r\ — 1. Then this 
assumption is 

(5.2) r = (n + r 2 )/2 > 1, r 2 - n = 1. 

With this condition, the Laplacian operator can be approximately expressed as 

d 2 1 d 2 



(5.3) 



A = 



dr 2 



de 2 ' 



With the boundary condition (5.1l and the operator (5.3 1, the eigenvalues and 
eignfunctions of the linear operator L\ = —A + B\ defined by (3.2 1 are given by 

p K (x) = K 2 (\~K 2 ), 

e K = cos kiOcos k 2 n(r — ri), 

e K — sin k\9 cos k 2 it(r — r%), 



K 



k\ 



ki,k 2 G Z. 
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Figure 5.1. All points in E A 
are singular points. 




Figure 5.2. The point A* is 
a singularity separation point 
from where two invariant sets 
T,\ and T\ are separated with 
S A being an attractor, and 
T A — > as A — > Aq. 



By (5.2 1, the first eigenvalue of L\ is 



1 



1 



/MA) = -2 (A -^), 
'0 '0 

which has multiplicity 2, and the first eigenfunctions are e\ — cos and e\ = sin ( 
Let Q, = S 1 x (ri,T2) satisfy 



Theorem 5.1 

hold true: 

(1) // 



5.2). Then the following assertions 



(2) 



73 > 



2rl 



then the phase transiti on o f (3.1 ) wi th (5.1) at Ao = 1/Vq js Type-I. Fur- 
thermore, the problem (3.1) with (5.1) bifurcates on A > Ao = r^ 2 to a cycle 
attractor = S* 1 which consists of singular points, as shown in Figure 5. 1 
and the singular points in can be expressed as 

1/2 



"A 



-1/2 



1 



A 2 cos(# 



A-i 



1/2 



1 



r73 



2 2 



4 ,u 6 

where Oq is the angle of u\ in Yi\. 
If 



73 < 



2rl 



"72: 



</ien i/ie transition is Type-II. Moreover, the problem bifurcates on A < Ao 
to a cycZe invariant set T\ = S 1 consisting of singular points, and there is 
a singularity separation at A* < Ao, where T\ and an S 1 attractor Ea = S 1 
are generated such that the system undergoes a transition at A = Ao from 
u = to u\ 6 as shown in Figure 
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Proof. Let v — ycosO + zsm8, u = v + $(y,z), and $ is the center manifold 
function. Then the reduced equations of ( |3.1| with (5.1 ) are given by 

dti 1 f ri+1 f 27T 

± = h.{\)y + - / / b 2 Av 2 + l3 Av :i + l2 Au 2 }e 1 1 drd9, 
at 7r j ri Jo 

dz 1 f Tl+1 f 2 ™ 

— = f3i(\)z + - / / [72AW 2 + 73 At) 3 + l2 &u 2 ]e\drdO. 

dt TT J ri J Q 



(5.4) 



Direct computation shows that (5.4 1 can be rewritten as 

c 2tt 



(5.5) 



dt TTTg 



dz 
dt 



Pi z 



1 



3tt 

T 

3tt 

T 



/■Z7T 

y(?/ 2 + z 2 ) + / 272$[j/cos 2 + zsin0cos0]e!0 
Jo 

7 3 z(y 2 + z 2 )+ / 272 $ [z sin 2 + ?/ cos sin 0]d0 
Jo 



and the center manifold function <I> = $(y, z) is 



272 



/3 2 (A)r 2 



(y 2 - z 2 ) cos 20 + 



472 



/3 2 (A)r 2 



yz sin 20, 



where 



W) = 4(A-^). 
'0 '0 



Putting $ into (|5.5|, we obtain the approximate equations of (5.4 1 as follows 



dy 



(5.6) 



dt 
dz 
dt 



1 



r 5 \4 



73 



2(A-4r - 2 ), 

s <- ( t + 7a ,.-ax > )*(y a + * 2 )- 



I /3 

/n \4 ' 2(A-4rn 2 ) 



At A = A = r 2 , we have 
3 

X 



7l 



2(A-4r " 2 ) 



3 -^22 
173- gr o72 . 



In the same fashion as used in Theorem 4.2 we derive from (5.6 1 the assertions 



of this theorem. Here the statement that the bifurcated cycle £>, = S 1 consists of 



singular points can proved by that the equation (3.1 1 with (5.1 1 is invariant under 
the translation 

tt(0,r)->u(0 + o ,r) V0 O e K 1 , 



which ensures that the singular points of (3.1 ) with (5.1 1 arise as a cycle S 1 . Thus, 
the theorem is proved. □ 



Now, we consider the problem that the equation (3.1) is defined in the whole 
space f2 = M. n (n > 2), with the periodic boundary condition 

(5.7) u(x + 2Kir) = u(x) WK = (fci, • • • , k n ) G Z n . 

In this case, the eigenvalues and eigenfunctions of L\ are given by 

/^(A) = |K| 2 (A-|X| 2 ), K=(k ir -- ,k n ), \K\ 2 = k 2 + --- + kl 



cos(fciiEi + • ■ ■ + k n x n ), e K = sin(fcixi + • • • + k n x n ). 
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It is clear that the first eigenvalue /3i(A) = A 
first eigenfunctions are 



1 of L\ has multiplicity 2n, and the 



ej = cos Xj , 



VI < j < n. 



Theorem 5.2. 



(1) If 



14 2 

73 > ^7 2 ; 



then the phase transition of (3.1) with (5.1) at Xq = 1 is Type-I. Moreover. 



(a) the problem bifurcates from (u, A) = (0, 1) to an attractor Y.\ homeo- 
morphic to a (2n — 1)- dimensional sphere 5 2ra_1 ; and 

(b) for each k (0 < k < n — 1), i/ie attractor T,\ contains C„- dimensional 
tori f n ~ k consisting of singular points. 

(2) // 

14 2 

73 < ^72 > 

i/ien the transition is Type-II. 
Proof. We only have to prove Assertion (2), as the remaining part of the theorem 



is essentially the same as the proof for Theorem 4.2 



Since the space of all even functions is an invariant space of L\ + G defined by 
(3.2 1, the problem (3.1) with (5.7) has solutions given in (4.31) with m = n in the 



space of even functions. 



1) and a 



By the translation invariance of (3.1 1 and (5.7), for each k (0 < k < n 
fixed index (Ji , 

(n — /c)-dimensional torus T n ~ k which consists of steady state solutions of (3.1 1 and 



(5.7 1. For example if (ji 



jk), the steady state solution associated with (4.31 1 gen erat es an 
~ k which ci 

,3k) = ' 



, k), the (n — fc)-dimensional singularity 



torus TP 



is 



T 



u(x + 1 



E 

j=k+i 



Vj cosOj + 0j) + o(\y\), V(0 r , 



where u(x) is the steady state solution of (3.1l with ( 5.7 1 associated with (4.31) 
with yi = • • • = y k = 0, y k+1 
Obviously, for a fix ed (J\ 



Vn 



■,3k) 



the 2 n ~ 



k steady state solutions of (3.1l and 
( |5.7| associated with ( |4.31| are in the same singular torus T™~ fe . 

,ik), 



Furthermore, 

for two different index fc-tuples (ji,-- - ,jk) and , i k ), the two associated 

singularity tori are different. Hence, for each < k < n — 1, there are exactly C k 
(n — fc)-dimcnsional singularity tori in £>,. Thus the proof is complete. □ 



6. Cahn-Hilliard Equations Coupled with Entropy 

When a phase separation takes place in a binary system, the entropy varies, and 
if the phase transition is Type-II, it will yield latent heat. Hence, it is necessary to 
discuss the equations (2.3 1, which are called the Cahn-Hilliard equations coupled 
with entropy. 
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To make the equations (2.3 1 non-dimensional, let 
x = lx' 

1 



k 2 



k 



-1 ,2 

/i = — I , 

k 2 



' ;4 
di = —I Oi, 



u = U U , 

I 2 



I 2 



k-2 
;2 



71 



-Sob , 



72 



-u b 2 



a 2 



73 



k 2 S 
I 2 



u a,2, 



S — SqS , 

I 2 

A = -H>i, 



-m 6 3 . 



Omitting the primes, equations (2.3 1 are in the following form 
dt 

at 



(6.1) 



= /iAS* — aiS 1 — a 2 w 2 , 
= -A 2 it- AAm + A(7i5u + 7 2 M 2 + 7 3 u a ), 
u{x, t)dx = 0, 



(9 



(u,Au,5)=0 on 90, 



<9n 

u(x, 0) = <p(x) 



By assumptions (2.2 1 and ( |2.4| , the coefficients satisfy 

^ > 0, ai > 0, a 2 > 0, 7i > 0, 73 > 0. 

Theorem 6.1. Let ft = n£ =1 (0,£ fc ) C R™ satwf/j/ that L = L x 
Lj(l < m < n) for any j > m. 

(1) For the case where m = 1, fe£ 



:73 



a 2 7i 



a 2 7l-£/ 



^ 72 

ai 2(aiL 2 + 4tt 2 ^) 3tt 2 



(a) If a < 0, then the phase transition of (6.1) at A = ir 2 /L 2 is Type-II 
and Assertion (1) in Theorem 4-1 holds true. 

(b) If a > the phase transition is Type-I and Assertion (2) in Theorem 
\4-l\ holds true. 

(2) For the case where m > 2, let 



9 



-73 - a 2 7i 



21? 



k £*i 2(a 1 L 2 + 4-7rV) a 1 L 2 + 27r 2 /^ 

(a) If a > 0, the phase transition of , 
Assertions (1) and (2) in Theorem 



13£ 2 7 2 2 
3tt 2 



6.1) at X = ir 2 /L 2 is Type-I and 
hold true. 



(b) If a < 0, i/ien i/ie phase transition is Type-II and Assertion (3) in 
Theorem\4-.2\ holds true. 



Proof. It suffices to compute the reduced equations of (6.1 ) on the center manifold. 
Similar to (4.181, the second order approximation of the reduced equation can be 
expressed as 

(6.2) 



dyi 

dt 



2J? 



[o-iVi +o-2Vi^2vj] 



2n 2 ™ 
Ll ...L n L^ y i 



<&i(y) cos — — cos —^-dx, 



L 



L 
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where (3\,a\ and oi are as in (4.18), and the center manifold function <&i(y) derived 
from the first equation in (6.1 ) can be expressed as 



(6.3) 



\K\>0 



®k(v) 



-OL 2 



AkII^II 2 



E 



Hi cos ■ 



L 



ipKdx 



where Xk and tpx are the eigenvalues and eigenfunctions of the following equation 

- [J,Aip K + OtlVK = >^kVk, 



dn 



-|an = 0, 



which are given by 



99 = 1, Vk = cos 
K=(k l /L 1 ,- 
Let 

Then we find 



kiirxi k n irx r , 
cos 



71 

k n /L n ), \Kf = J2k 2 /L 2 



for fcj 6 Z, 1 < i < n. 



i^i = (6a/ Li, • • • , 6 in /L n ). 



$0 =- 



-«2 



aiLi ■ • • L r , 



E 



J/j cos 



L 



dx = 



Ei 



-CK2 



2ai 



-a 2 

akIi^ii 2 y n 





E 



J/iCOS 



rf.r 



if K^Ki + K r , 



~a 2 



2 E f'fj / cos cos 1 —r L PKdx a k = K t + K r , 



hJ —i 

for some 1 < r, I < m. Thus we derive that 

-a 2 



A" 



AaII^aII 2 ^ 

XkW^kW 2 ^ Jr: 
-Q.1 



!i~ ! cos 2 ^efc 



2(ai + 4tt 2 ^/ j L 2 ) 
-2a 2 



2tt 2 h/L 2 



HK = K j +K j , l<j<m 

I cos — - — cos —j^-ex u K = Ki + Kj, i f= j 
Jn L L 

ifK = 2Kj, l<j<m, 
if K = Ki + Kj , ij^j, 1 < i, j < m. 
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Putting <j> and <&k in fl6.3| ) we obtain 

n n _ 

*i =^7 L »j " 2(ai + 47r W l 2 ) 2. »j cos — 

3=1 ' ' 3=1 

2a 2 ""^i 7ra; 7 

" ai + 27T ^ ^ "i" ^' 

i<j 

Then, inserting <I>i into (|6.2|), we derive the following reduced equations: 

_ QgTi Q27i \ 3 

CTl ax 2(ai +4vr 2 /i/i 2 y tfi 



(6.4) f=A(Ate,-^ 



a 2 7x 2a 2 7x 



(r2 -V"a 1 + 2,V/i 2,!/i;j 



2 
3 



1 < i < m. 



Then the remaining part of the proof can be achieved in the same fashion as the 
proofs for Theorems |4.1| and |4.2| The proof is complete. □ 



Remark 6.1. From the phcnomcnological viewpoint, the coefficient fi > in (6.1 ) 



is small. If let /i = 0, then in the equilibrium state we have that S = —^u 



In this case, (6.1) are referred to the original Cahn-Hilliard equation (3.1l with 



73 = 73 — a27x/«i as the coefficient of the cubic term, and the criterion a = in 



Assertions (1) and (2) of Theorem 6.1 are respectively equivalent to 



, 2L^i , 26^72 

73=9^2" ^ 73=^^2", 

which coincide with these in Theorems |4.1| and |4.2| Hence, if we consider the 
coefficient /i > small, then the criterion a — are respectively equivalent to 

, _ 2 L 2 7| 4 7r 2 ^q 2 7i , _ 26 L 2 j 2 4 7r 2 /za 2 7i 

73 ~ 9 ir 2 3 a 2 L 2 ' 73 "27 tt 2 3 a\I? ' 

Hence the item — (47r 2 /ia27i)/3af L 2 is the effect yielded by /j,AS. □ 

7. Physical remarks 

We now address the physical significance for the phase transition theorems ob- 
tained in the previous sections. 

7.1. Equation of critical parameters. For a binary system, the equation de- 
scribing the control parameters T,p, f2 at the critical states is simple. 

We first consider the critical temperature T c . There are two different critical 



temperatures T\ and T in the Cahn-Hilliard equation. Ti is the one given by (2.8 I 
at which the coefficient bi(T,p) or A — —l 2 bi(T,p)/k will change its sign, and T( 
satisfies that Ao < Ai, and for fixed p, 

<pi ifT>T , 

(7.1) ACT] { = P i if T = T , 

>pi ifT<T , 



where pi is the first eigenvalue of (3.4), which depends on the geometrical properties 
of the material such as the size of the container of the sample f2. When £1 — 
(0,L) m x D C K." is a rectangular domain with L > diameter of D,pi = tt 2 /L 2 . 
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Hence, in general at the critical temperature T\ a binary system does not undergo 
any phase transitions, but the phase transition does occur at T = T . 
At T\ and Tq we know that 

(7.2) A(Ti)=0, A(T )=pi. 

For a rectangular domain, pi = ir 2 /L 2 , therefore from (7.2 1 we see that T\ is a limit 



of the critical temperature To of phase transition as the size of tends to infinite. 

In fact, for a general domain, it is easy to see that the first eigenvalue pi of the 
Laplace operator is inversely proportional to the square of the maximum diameter 

of n-. 

(7-3) pi ~ 

where L represents the diameter scaling of £1. 



Thus the equation of critical parameters in the Cahn-Hilliard equation, by (7.2) 
and ( |7.3| , is given by 

(7.4) A(T,p) = -^, 

where C > is a constant depending on the geometry of J7. According to the 
Hildebrand theory (see Reichl 7 ) , the function X(T, p) can be expressed in a explicit 
formula. If regardless of the term |Vm| 2 , the molar Gibbs free energy takes the 
following form 

(7.5) g — pa(1 — u) + pbu + RT{\ — u) ln(l — u) + RTuliiu + au(\ — u), 

where pa , P b are the chemical potential of A and B respectively, R the molar gas 
constant, a > the measure of repel action between A and B. Therefore, the 



coefficient b\ in (2.7 1 with constant p is 



h = 7rf|«=u = — 77- zRT-a (a = a(p)), 

au z uo(l — uo) 

where u — ug is the constant concentration of B. Hence 

w , l\ 2al 2 l 2 R 

\(T,p) = --b 1 = — - —— -T. 

k k fcuo(l — uo) 



Thus, equation (7.4) is expressed as 



(7 6) — T 

kuo(l — Uo) K L 2 



Equation (7.6) gives the critical parameter curve of a binary system with constant 



pressure for temperature T and diameter scaling L of container fl. Because T > 0, 



from (7.5) we can deduce the following physical conclusion. 



Physical Conclusion 7.1. Under constant pressure, for any binary system with 
given geometrical shape of the container f2 ; there is a value L > such that as the 
diameter scaling L < Lq, no phase separation takes place at all temperature T > 0. 
and as L > Lq phase separation will occur at some critical temperature Tq > 



satisfying (7.6) 



We shall see later that it is a universal property that the dynamical properties 
of phase transitions depend on the geometrical shape and size of the container or 
sample il. 
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7.2. Physical explanations of phase transition theorems. We first briefly 
recall the classical thermodynamic theory for a binary system. Physically phase 
separation processes taking place in an unstable state are called spinodal decompo- 
sitions; see Cahn and Hilliard pQ and Onuki [B]. When consider the concentration 
u as homogeneous in f2, then by (7.5 1 the dynamic equation of a binary system is 
an ordinary differential equation: 

du dg 
~dt ~ 



(7.7) 



du 



2au- i?Tln 



1 - u 



VA- IJ-B ~ a. 



Let Mo (0 < i*o < 1) be the steady state solution of (7.7 1. Then, by the Taylor 
expansion at u = no, omitting the nth order terms with n > 4, (7.7 1 can be 
rewritten as 

- Xv + b 2 v 2 + b 3 v 6 , 



(7.8) 
where 



dt 



V = U — Uq , 

1 - 2u 



A = 2a - 



u (l - u ) 



RT, 



2ui(\-u y 

It is easy to see that 



RT, 



1 1 - uq - 2u\ + 3ul 
3 ug(l-Mo) 4 



RT. 



= if u = 

^0 if^o^, 

b 3 < V0 < u < 1. 

It is clear that the critical parameter curve A = in the T — uq plane is given by 

T = 2au (l - u )/R, 

which is schematically illustrated in the classical phase diagram; see the dotted line 
in Figure 7.1 We obtain from (7.8) the following transition steady states: 

+ -1 



26a 



(6 2 ± \ b\ - 4fe 3 A) for b\ - 46 3 A > 0. 



By Theorem 
namely 



3.2 



(7.9) 



T* = 



we see that there is T* = T*(u Q ) satisfying that b\ - 46 3 A = 0; 
2au (l-Uo) a , 3(1 - 2w ) 2 (l - uo) 



R(l-f3(u )) : 



0{u o ) 



16(1 - u - 2u\ + 3u%) ' 



such that if T n < T < T* 



v + is 



locally stable (metastable) 
unstable 

locally stable (metastable) 
unstable 



for < u < 
for - < u < 1, 
for i < uq < 1, 
for < u < 
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stable region 



Figure 7.1. Typical 
phase diagram from 
classical thermody- 
namic theory. 




Figure 7.2. The state u = 
M B is stable if T c = T < T, 
and the state uo is unsta ble i f 
T < T , where T is as in frill. 




Here 



T K) = - > ToOo) 

1 - Po{u ) 



is illustrated by the solid line in Figure 7.1 This shows that the region Tq(uq) < 
T < T*(uq) is metastable, which is marked by the shadowed region in Figure |7T] 
See, among others, Reichl [7], Novick-Cohen and Segal [5] , and Langer [2] for the 
phase transition diagram from the classical thermodynamic theory. 

In the following we shall discuss the spinodal decomposition in a unified fashion 
by applying the phase transition theorems presented in the previous sections. 

As mentioned in the Introduction, phase separation processes of binary systems 
occur in two ways, one of which proceeds continuously depending on T, and the 
other one does not. Obviously, the classical theory does not explain these phenom- 
ena. In fact, the first one can be described by the Type-I phase transition, and the 
second one can be explained by the Type-II and Type-Ill phase transitions. 

We first consider the case where the container f2 = n" =1 (0,Lj) with L = L\ = 



4.2 



• • • = L m > Lj(j > to) is a rectangular domain. Thus, by Theorems 4.1 and 
(or Theorem |6.1[ ) there are only two phase transition types: Type-I and Type-II 
with the type of transition depending on L. We see that if 

9 7T 2 7 3 



2 it 
27 7T 2 7 3 



for m = 1, 
for to > 2, 



26 7 | 

then the transition is Type-I, i.e., the phase pattern formation gradually varies as 
the temperature decreases. In this case, no meta-stable states and no latent heat 



appear. The phase diagram is given by Figure 7.2 where the solid lines uj (i = 1, 2) 
represent the transition solutions. 
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u 



Figure 7.3. A Type-II phase transition. 



If L satisfies that 



L 2 > 



9 tt 2 7 3 
2 ll 
27^73 
26 7 | 



for m = 1, 



for m > 2, 



then the phase transition is Type-II. Namely, there is a leaping change in phase 
pattern formation at the critical temperature T c . The phase diagram for Type-II 



transition is given by Figure 7.3 



In Figure 7.3 To is the critical temperature as in (7.1 I, T* is defined by (7.9) and 
is the saddle-node bifurcation point of (7.8). The constant concentration u = u is 



stable in T c < T, is meta-stable in T < T < T c , and is unstable in T < T . The 
two bifurcated states Uf and U% from T* are meta-stable in Tq < T < T* , and 



are stable in T < T . Here for i = 1, 2, Uj 



Mo, and Ui are the separated 



solutions of (3.1 1 with (3.5 1 from T 



There is a remarkable difference between Type-I and Type-II transitions. The 
Type-I phase transition occurs at T — T and Type-II does in T < T < T*. 
Furthermore, latent heat is accompanied the Type-II phase transition. Actually, 
when a binary system undergoes a transition from uo to Uf [i — 1,2), there is a 
gap \Uf - u \ 2 = \uf\ 2 > e > for any T < T < T* . By the first equation in Kty 



it yields a jump of entropy between uq and Ui 



SSi = / Sdx — 



OL 2 

ax 



T\2 



dx < 0, 



where S = Si 
is given by 



in "i Jn 

Sq represents the entropy density deviation. Hence the latent heat 



SH = T8S l 



a 2 T 



,T|2 



TT 



dx < 0, 



which implies that the process from uq to U( is exothermic, and the process from 
Uf to uo is endothcrmic. 

Now, we consider the case where the container is non rectangular. Thus, by 
Theorem |3.2| the transition is Type-Ill, and its phase diagram is given by Figure 
PI 



In Figure 7.4 To and T* are the same as those in Figure 7.3 The state uq is 
stable in T* < T, is metastable in T < T < T* , and unstable in T < T . The 
equilibrium state U± separated from T* is metastable in To < T < T* , and is stable 
in T < T . However, the equilibrium state U% separated from T* is unstable in 
T < T < T*, and is metastable in T < T . 
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Similar to the Type-II, the Type-Ill phase transition has also latent heat, which 
occurs in T < T < T*. But the difference between Type-II and Type-Ill is that 
Type-II has 2m (to > 1) stable equilibrium states separated from T — T* , but 
Type-Ill has just one. The 2m stable states of a Type-II transition are of some 
symmetry caused by f2, and we shall investigate it later. A particular aspect of 
Type-Ill is that there is a state bifurcated from (u,T) = (uo,Tq), which is 
rarely observed in experiments. 



7.3. Symmetry and periodic structure. Physical experiments have shown that 
in pattern formation via phase separation, periodic or semi-periodic structure ap- 
pears. From Theorems |5.1| and |5.2| we see that for the loop domains and bulk 
domains which can be considered as M. n or K m x D (D C K n_m ) the steady state 
solutions of the Cahn-Hilliard equation are periodic, and for rectangular domains 
they are semi-periodic, and the periodicity is associated with the mirror image 
symmetry. 

Let n = (0,L) m x D (to > 1). By Remark 
Actually, in this case the following space 



4.1 



Theorem 



4.1 



is valid for f2. 



H = < lie L 2 (n)\ u = 



oo 7 

y k cos 

\K\ = 1 



■ cos 



L 



C H 



is invariant for the Cah n-Hilliard equation (3.1 1 and (3.5 1. All separated equilibrium 
states in Theorem ■ 



4.1 



are in H . From the physical viewpoint, all equilibrium states 
u{x) and their mirror image states u(L—x', x") are the same to describe the pattern 
formation. Mathematically, under the mirror image transformation 

X ► (Z/ X , X ) , X — (xi , • • • , X m ) , X (l£m+l i ' ' ' i ^n) ? 

the Cahn-Hilliard equation ( |3.1[ ) with (3.5 1 is invariant. Hence, the steady state 
solutions will appear in pairs. In particular, for Type-I phase transition, there is a 
remarkable mirror image symmetric. We address this problem as follows. 

Let m = 1 in ft = (0,L) m x D. By (4.2 1 there are two bifurcated stable equilib- 
rium states, and their projections on the first eigenspace are 



U\ = y cos 



TTXl 

17 



u 2 



-ycos 



TTXl 



It is dear that u 2 (xi) 



y= v/2(A-^/^)/(J 73 

= m(L - xi). 



3tt2 
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Let m = 2. By Theorem |4.2| the bifurcated attractor contains 8 equilibrium 
states, whose projections are given by 



Ui =±J/oei, u 2 = ±yi(ei + e 2 ), 
u f = ±yoC2, uf = ±yx(-e 1 + e 2 ), 

where ei = (cos and e 2 = (0,cos form an orthogonal basis in IR 2 , and 

?/i = \/2(X — 7r 2 / 'L 2 ) I '\/<7\ + (72 with <7i,<72 as in (4.191. 

The se eight equilibrium states constitute an octagon in W 2 , as shown in Figure 
7.5 and they are divided into two classes: A\ = {uf,uf } and A3 = {w 2 , u 4 } 
by the ^-rotation group G(f ). Namely, with the action of G(^), A (i = 1,2) are 
invariant: 

BueA l: Vu G A and £ e <?(-), 
where G(|) consists of the orthogonal matrices 

*-*G !)■ "MS "M! ^ 

The stability of the equilibrium states Uk is associated with both classes Ai and 
A2, i.e., either the elements in Ai are stable, or those in A2 are stable; see Figure 



4.2 By (4.331 we can derive the criterion as follows 

± 22 L ^2 2 6 i 2 7l . , 

u 2k+i e -^i are stable <=> ^ > 73 > ™ — for fc = 0, 1, 

9 7T 27 7T 

+ , 22 L 2 7 2 
uf k e A are stable 73 > — ^— for k — 1,2, 



In Figure 7.5 we see that elements in Ai and A have a J difference in their 
phase angle. However, in their pattern structure, u^ K+1 £ A and uf k £ Ai also 
have a 7r/4 deference at the angle between the lines of u 2K+1 = and uf K = 0. In 
fact, the lines that 

uj = —U-y = Uq cos — — — U, and 

u 3 = ~ u 3 = VO COS — — = 



32 



MA AND WANG 



are given by x\ = L/2 and X2 = L/2 respectively, as shown in Figure [7^ a), and 
the lines 



1*2 = — u 2 = yi I cos 



TTXl 



cos ■ 



TTX 2 



= 0, and 



L L 

x 4 — -u 4 = yi y— cos — - + cos —j— J = 

are given by X2 = L — x\ and X2 — x\ respectively as shown in Figure |7l3[b). 



L /2 

(a) 




Figure 7.6. 



Let m = 3. Then the bifurcated attractor contains 26 equilibrium states 
which can by divide into three classes by the 3-dimensional (§, §, f ) -rotation group 
G (§, f , f) as follows 

Ai = {uf = ±y ex,uf = ±y e 2 ,u^ = ±y e 3 } , 

Ai = {uj = ±yi (ei + e 2 ),w* = ±t/i (e x + e 3 ),itg = ±J/i(e 2 + e 3 ), 

"7 = ±2/1 (~ei + e 2 ),u* = ±2/1 (-e x + e 3 ),itg = ±2/1 (-e 2 + e 3 )}, 
■^3 = {w^q = ±2/1 (ei + e 2 + e^it^ = ±j/i(-ei + e 2 + 3 3 ), 

u i2 = ±2/i ( e i - e 2 +e 3 ),itf 3 = ±2/1 (ei + e 2 - e 3 )}. 

Only these elements in Ai or in A3 are stable, and they are determined by the 
following criterion 



22 i 2 7 2 



elements in Ai is stable <^> 



26 T 2 7 2 



> 73 > 



elements in A3 is stable 7 3 > 



9 7T 2 M 27 7T 2 

22 L 2 7 | 



9 rr 2 



7.4. Critical exponents. From (4.2 1 and (4.341 we see that for Type-I phase 
transition of a binary system the critical exponent /3 = h. In this case, it is a 
second order phase transition with the Ehrcnfcst classification scheme, and there is 



a gap in heat capacity at critical temperature To- To see this, by (4.2 1 and (4.34) 
we have 

f if T < T, 

U ~ I a(A(T)-7r 2 /L 2 ) 1 / 2 e 1 +o(|A-rr 2 /L 2 | 1 / 2 ) if T < T , 
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and the free energy for (3.1) at u T is 
1, 



|v*T-£l«T + o(l«T) 



dx 



A« r ]u T +^(u r ) 3 +i7!i(« T ) 4 ^ 



if T < T, 



2^2 



-(A-TrVi 2 ) 



[e? + ^73^ + o(|A - 7r 2 /i 2 | 2 ) if T > T . 



Thus, the heat capacity C at T = Tq satisfies 



C+ -C~ 



-T 



d 2 F{u T ) 



dT 2 



T 



d 2 F{u T ) 



dT 2 



aiT 



dXY 
dTj 



T=T n 



It is known that dX/dT ^ 0; hence the heat capacity at T = T has a finite jump. 



From (3.141 we know that for the Type-Ill case, the critical exponent (3 = 1. 



Thus, it is not hard to deduce that the continuous phase transition in Typc-III is 
of the 3rd order. 

Appendix A. Dynamic Transition Theory for Nonlinear Systems 

In this appendix we recall some basic elements of the dynamic transition theory 
developed by the authors [3JH], which are used to carry out the dynamic transition 
analysis for the binary systems in this article. 

A.l. New classification scheme. Let X and X\ be two Banach spaces, and 
X\ C X a compact and dense inclusion. In this chapter, we always consider the 
following nonlinear evolution equations 

(A.l) Tt =L x u+G{u,\), 

u(0) = ip, 

where u : [0,oo) — > X is unknown function, and A £ K 1 is the system parameter. 

Assume that L\ : X\ — > X is a parameterized linear completely continuous field 
depending continuously on A e I 1 , which satisfies 

L\ = —A + B\ a sectorial operator, 
(A. 2) A : X\ — » X a linear homeomorphism, 

B\ : X\ — > X a linear compact operator. 

In this case, we can define the fractional order spaces X a for a £ K 1 . Then we also 
assume that £?(•, A) : X a — > X is C r (r > 1) bounded mapping for some < a < 1, 
depending continuously on A € M 1 , and 

(A.3) G(u,A) = o(||u|| x J, VAeR 1 . 



Hereafter we always assume the conditions (A. 2 1 and (A.3 1, which represent that 



the system (A.l I has a dissipative structure. 



Definition A.l. We say that the system {A.l) has a transition of equilibrium from 
(it, A) = (0, Aq) on A > Aq (or A < Aq) if the following two conditions are satisfied: 



(1) when A < A (or A > \q),u = is locally asymptotically stable for {A.l) 
and 
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(2) when A > Ao (or A < Ao), there exists a neighborhood U C X of u = 
independent of A, such that for any if S U \ T\ the solution U\(t,ip) of 
(A.l) satisfies that 

limsup ||u A (*,¥>)IU > S(X) > 0, 

lim 5(A) > 0, 

A^A 

where T\ is the stable manifold of u = 0, loii/i codim T\ > I in X for 
A > A for A < A ). 



Obviously, the attractor bifurcation of (A.l) is a type of transition. However, 
bifurcation and transition are two different, but related concepts. Definition |A.1| 
defines the transition of (A.l I from a stable equilibrium point to other states (not 
necessary equilibrium state) . In general, we can define transitions from one attrac- 
tor to another as follows. 

Let the eigenvalues (counting multiplicity) of L\ be given by 



Assume that 



(A.4) 



(A.5) 



Re A (A) 



Re (3j(\ ) < 




{/?i(A)eC | j \. ■>.■■■} 



if A < Ao, 
if A = A , 
if A > A , 



VI < i < m, 



Vj > m + 1. 



The following theorem is a basic principle of transitions from equilibrium states, 
which provides sufficient conditions and a basic classification for transitions of non- 
linear dissipative systems. This theorem is a direct consequence of the center man- 
ifold theorems and the stable manifold theorems; we omit the proof. 



Theorem A.l. Let the conditions {A. 4) and (A.5) hold true. Then, the system 
(A.l) must have a transition from (u, A) = (0, Ao), and there is a neighborhood 
U C X of u = such that the transition is one of the following three types: 

(1) CONTINUOUS Transition: there exists an open and dense set U\ C U 
such that for any if G U\, the solution U\(t, if) of (A.l) satisfies 

lim limsup ||wa(£, <p)\\x = 0. 

A — >Ap t — >oo 



In particular, the attractor bifurcation of (A.l) at (0, Ao) is a continuous 
transition. 

(2) Jump Transition: for any Ao < A < Ao + e with some e > 0, there is an 
open and dense set U\ C U such that for any if G U\, 

limsup 1 1 Ua (*,</>) | U > 5 > 0, 

t — >oo 

where S > is independent of X. This type of transition is also called the 
discontinuous transition. 

(3) Mixed Transition: for any A < A < A + e with some e > 0, U can be 
decomposed into two open sets and f/^ 1 (Uf not necessarily connected): 



vl n c/ 2 A = 
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such that 



lim limsup \\u(t, <p)\\x = C 
limsup </?)||y > <5 > 

t — >oo 



V<y9 g t/ 2 A 



With this theorem in our disposal, we are in position to give a new dynamic 
classification scheme for dynamic phase transitions. 

Definition A.l (Dynamic Classification of Phase Transition). The phase transi- 
tions for {A.l) at X = Xq is classified using their dynamic properties: continuous, 
jump, and mixed as given in Theorem \ A.1\ which are called Type-I, Type-II and 
Type-Ill respectively. 

An important aspect of the transition theory is to determine which of the three 
types of transitions given by Theorem |A. l| occurs in a specific problem. By reduction 
to the center manifold of (A.l), we know that the type of transitions for (A.l I at 
(0, Ao) is completely dictated by its reduction equation near A = Ao, which can be 
expressed as: 

dx 



(A.6) 



dt 



= J\x + PG(x + $(x, A), A) for x G 



where J\ is the m x m order Jordan matrix corresponding to the eigenvalues given 
by (A. 4), Q(x, A) is the center manifold function of (A.l I near Ao, P : X — > E\ is 
the canonical projection, and 

E x = Ui<i< m U feeJV {u G X\ | (L x - fr(X)) k u = 0} 

is the eigenspace of Lx ■ 

By the spectral theorem, (A.6 1 can be expressed into the following explicit form 

(A.7) 



dx 



Jxx + g(x, A), 



where 
(A.8) 



g{x,X) = (siO,A),-- - ,g m (x, A)), 

m 

9j{x, X) =< G(^2x t e. t + $(x,A),A),e* > VI < j < m. 



Here and e* (1 < j ' < m) are the eigenvectors of Lx and L* x respectively corre- 
sponding to the eigenvalues (3j(X) as in (A.4.I. 

In particular, if G{u, X) has the Taylor expansion 

(A.9) G(u,X) = G k (u,X)+o(\\u\\ k x ), 



for some k > 2, where Gk(u, X) is a k- multilinear operator, then (A.7 1 can be 
rewritten as 

dx 



(A.10) 
where 



dt 



= Jxx + 9k(x 1 X) + o(\x\ k ), 



gic(x,X) = (g kl (x,X),- ■ ■ ,g km (x,X)), 

n 

g kj (x,X) =< G k (^2xiei,X),e* > 



VI < j < m. 



i=l 



3(5 



MA AND WANG 



When x = is an isolated singular point of gk(x, A), in general the transition of 



(|A.1|) is determined by the first-order approximate bifurcation equation of (A. 10) 

dx 



as follows 
(A.ll) 



dt 



= J\x + g k (x, A). 



The following theorem is useful to distinguish the transition types of (A.ll at 
(«,A) = (0,Ao). 



Theorem A. 2. Let the conditions (A. 4) and (A. 5) hold true, and U C R m be a 
neighborhood of x — 0. Then we have the following assertions: 



(1) If the transition of (A.l) at (0,Ao) is continuous, then there is an o pen a nd 
dense set U C U such that for any x E U the solution x(t,Xo) of (A.l) at 
A = Ao with x(0,xq) = xq satisfies that 

lim x(t, Xq) = 0. 

t — KX) 

(2) If there exists an open and dense set U C U such that for any Xq E U the 
solution x(t,Xo) of (A.l) at X = Xq satisfies 

lim sup \x(t, xq) \ 7^ 0, 

t — >oo 

then the transition is a jump transition. 

(3) // the transition is mixed, then there e xists an open set U C U such that 
for any xq € U the solution x(t, xq) of (A.l) at X = Xq with x(0, Xq) = xq 
satisfies 

lim x(t, Xq) = 0. 

t — >oo 



(4) If the vector field in ![A.1\ ) at X = Xq satisfies 

< J\„x + g{x, Aq), x >< Vie G U, x ^ 0, 



then the transition of (A.l) at (0, Ao) is an S m 1 -attractor bifurcation, and 
the transition is continuous. 



(5) If the vector field g(x,Xo) given by (A. 8) satisfies 

<g(x,Xo),x»0 Vx G U, x ^ 0, 
then the transition is jump. 



In general, the conditions in Assertions (l)-(3) of Theorem A. 2 are not sufficient. 
They, however, do give sufficient conditions when (A.l I has a variational structure. 
To see this, let (A.ll be a gradient-type equation. Under the conditions (A. 4 1 and 
(A. 5 1, in a neighborhood U C X of u = 0, the center manifold M c in U at A = A 
consists three subsets 

M c = w u + W s + D, 

where W s is the stable set, W u is the unstable set, and D is the hyperbolic set of 
(A. 7 1. Then we have the following theorem. 



Theorem A. 3. Let (A.l) be a gradient-type equation, and the conditions (A. 4-) 
and (A. 5) hold true. If u = is an isolated singular point of (A.l) at X = Xq, then 
we have the following assertions: 
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(1) The transition of (A.l ) at (u, A) = (0, Ao) is continuous if and only if u = 
is locally asymptotically stable at X — X , i.e., the center manifold is stable: 
W s . Moreover, (A.l) bifurcates from (0, Ao) to minimal attractors 



M c 

consisting of singular points of (A.l) 



(2) // the stable set W s of (A.l) has no interior points in M c , i.e., M c 



W u + D, then the transition is jump. 



A. 2. Transitions from simple eigenvalues. We consider the transition of (A.l ) 



from a simple critical eigenvalue. Let the eigenvalues f3j(X) of L\ satisfy (A.4) and 
( |A.5[ ) with m = 1. Then the first eigenvalue /3i(A) must be a real eigenvalue. Let 
ei(A) and (A) be the eigenvectors of L\ and L* x respectively corresponding to 
01(A) with 



L Ao ei=0, L* x el=0, <e 1 ,e* 1 >=l. 



Let $(a;, A) be the center manifold function of (A.l I near A = A . We assume that 



(A.12) 



< G(xex + $(sg, A ), A ), e* >= ax k + o(\x\ k ), 



where k > 2 an integer and a/Oa real number. 













V/- 


u = 





(a) 



(b) 



Figure A.l. Topological structure of the jump transition of (A.l I 
when k=odd and a > 0: (a) A < A ; (b) A > A . Here the 
horizontal line represents the center manifold. 



(a) 



Figure A. 2. Topological structure of the continuous transition of 
(|A.1| when /c=odd and a < 0: (a) A < A ; (b) A > A . 
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(a) 



(b) 



(c) 



Figure A. 3. Topological structure of the mixing transition of 
( A. 1 1 when fc=even and a ^ 0: (a) A < Ao; (b) A = A ; (c) 
A > Aq. Here U* is the unstable domain, and 17% the stable do- 



Theorem A. 4. Assume (A. 4.) and (A. 5) with m = 1, and (A. 12). If k — odd and 



a =/= in (A. 12) then the following assertions hold true: 



(1) If a > 0, then (A.l) has a jump transition from (0, Ao), and bifurcates on 
A < Ao to exactly two saddle points v* and v* with the Morse index one, 

1X71 



as shown in Figure 



(2) lfa<0, then (A.l) has a continuous transition from (0,Ao), which is an 
attractor bifurcation as shown in Figure \A.2\ 

(3) The bifurcated singular points v* and v* in the above cases can be expressed 
in the following form 

v* = ±|/3 1 (A)/a| 1 / fc - 1 ei (A) + odfrl 1 /*-*). 



Theorem A. 5. Assume (A.Jty and \A.5 ) with m = 1, and (A. 12). If k — even and 
a =/= 0, then we have the following assertions: 



(1) (A.l) has a mixed transition from (0,Ao). More precisely, there exists a 
neighborhood U C X of u = such that U is separated into two disjoint 
open sets IF* and t/^ by the stable manifold T\ of u = satisfying the 
following properties: 

(a) U = U* + U*+T X , 

(b) the transition in U* is jump, and 

(c) the transition in is continuous. The local transition structure is as 



shown in Figure A. 3 



(2) {A.l) bifurcates in 11% to a unique singular point v on A > Aq, which is 



an attractor such that for any ip G U : 



2 1 

lim \\u(t,tp) -v x \\x = 0, 

t — >oo 



wher e u(t,tp) is the solution of (A.l). 
(3) (A.l ) bifurcates on A < Aq to a unique saddle point v* with the Morse index 



one. 



(4) The bifurcated singular point v* can be expressed as 



v* = -(/5 1 (A)/«) 1 /( fc - 1 ) ei + odAlVC*-*)). 

A. 3. Singular Separation. In this section, we study an important problem asso- 



ciated with the discontinuous transition of (A.l), which we call the singular sepa- 
ration. 
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Definition A. 2. (1) An invariant set £ of (A.l) is called a singular element 
if E is either a singular point or a periodic orbit. 



(2) Let Si C X be a singular element of (A.l) and U C X a neighborhood of 
El. We say that (A.l) has a singular separation of 'E at A = Ai if 



(a) (A.l) has no singular elements in U as X < Ai (or X > X\), and 
generates a singular element Ei C U at A = Ai, and 

(b) there are branches of singular elements £>, which are separated from 
Ei for A > Ai (or X < Xi), i.e., 

lim max disttx, Ei) = 0. 

A special case of singular separation is the saddle-node bifurcation denned as 
follows. 



Definition A. 3. Let ui G X be a singular point of {A.l) at X = X\ with u\ ^ 0. 
We say that (A.l) has a saddle-node bifurcation at («i,Ai) if 

(1) the index of L\ -\- G at (ui,X\) is zero, i.e., ind(—(L\ 1 +G),ui)=0, 



(2) there are two branches T 'i(A) andT2(X) of singular points of (A.l), which 



Ai. 



are separated from u\ for X > Ai (or X < X±), i.e., for any u\ € Ti(A) 
(i = 1,2) we have 

u\ — > u\ in X as X 

and 

(3) the indices of u\ € Tj(A) are as follows 
ind{-{L x + G),u x ) = \ 1 



1 



ifu x e r 2 (A), 
ifuxer^x). 



Intuitively, the saddle-node bifurcation is schematically shown as in Figure |A.4| 
where the singular points in Ti(A) are saddle points and in ^(A) are nodes, and 
the singular separation of periodic orbits is as in shown Figure |A.5| 




Figure A. 4. Saddle-node bifurcation. 

For the singular separation we can give a general principle as follows, which 
provides a basis for singular separation theory. 



Theorem A. 6. Let the conditions (A. 4) and \ A.5) hold true. Then we have the 
following assertions. 

(1) If (A.l) bifurcates from (it, A) = (0, Aq) to a branch Ea of singular elements 



on X < X which is bounded in X x (-co, Ao) then \A. lty has a singular 
separation of singular elements at some (E ,Ai) C X x (— oo,Aq). 
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Figure A. 5. Singular separation of periodic orbits. 



(2) If the bifurcated branch T.\ consists of singular points which has index —1, 
i.e., 

ind(-(L x + G),u\) = -1 Vu A e E x , A < A , 

then the singular separation is a saddle-node bifurcation from some (u\, X\) G 
X x (—oo, Aq). 



H,Xi 

H, we assume that A: Hi - 



We consider the equation ( A.l I defined on the Hilbert spaces X 
Let L\ = —A + \B. For L\ and G(-,X) : H 1 
symmetric, and 

(A. 13) < Au,u>h>c\\u\\ 2 Hi/2 , 

(A. 14) < Bu,u > H > c\\u\\ 2 H , 

(A. 15) < Gu,u > H < -ci\\u\\ p H + c 2 \\u\\ 2 H , 

where p > 2, c, ci, C2 > are constants. 



H is 



Theorem A. 7. Assume the conditions (A. 3), (A. 4) and (A.13)-(A.15), then (A.l) 
has a transition at (it, A) = (0, Ao), and the following assertions hold true: 

(1) lfu = is an even-order nondegenerate singular point of L\ + G at A = Ao . 
then (A.l) has a singular separation of singular points at some (ui,Ai) € 



H x (— oo, A ). 



(2) If m = 1 and G satisfies (A. 12) with a > ifk=odd and a ^ if k=even, 
then (A.l) has a saddle-node bifurcation at some singular point (ui,Ai) 
with Ai < Ao- 

A. 4. Transition and Singular Separation of Perturbed Systems. We con- 
sider the following perturbed equation of ( |A.1[ ): 

du 



(A.16) 



dt 



= (L X + S s x )u + G(u,\)+T e (u,\), 



where L\ and G\ are as in (A.l), S\ : X a — » X is a linear perturbed operator, 
: X a — > X a C 1 nonlinear perturbed operator, and X a the fractional order 
space, < a < 1. Also assume that G\,T^ are C 3 on w, and 



(A.17) 



IISaII <e, 

T e {u, A) = o(||u||jcJ. 



Let (A.4| and (A. 5) with m = 1 hold true, G(w, A) 



G 2 (u,A)+o(||u||V), where 



G2(-, A) is a bilinear operator, and 

(A.18) b=< G 2 (e,X ),e* >^ 0, 
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where e£l and e* £ X* are the eigenvectors of L\ and L* x corresponding to /3i(A) 
at A = A respectively. 

We now consider the transition associated with the saddle-node bifurcation of 
the perturbed system (A. 16 1. Let h(x,X) be the center manifold function of (A.l) 
near A 



(A.19) 



Ao- Assume that 

< G(xe + h(x, X Q ), A ), e* >= hx 3 + o(\x\ 3 ), 



where bi ^ 0, and e and e* are as in (A. 18 1. 
Then we have the following theorems. 



Theorem A. 8. Let the conditions (A. 4) and (A. 5) with m = 1, and (A.19) hold 
true, and &i < 0. Then there is an e > such that if iS| and satisfy (A. 17), 
then the transition of (A. 16) is either continuous or mixed. If the transition is 
continuous, then Assertions (2) and (3) of Theorem A. 4 are valid for (A. 16). If 
the transition is mixed, then the following assertions hold true: 

(1) (A. 16) has a saddle-node bifurcation at some point {u\, \\) £Xx (-co, Ag), 
and there are exactly two branches 

r A = {(u A ,A)| Ai<A<Ag + 5} i = 1,2, 
separated from (iti,Ai) as shown in Figure \A~4~ which satisfy that 

V(u A ,A)er A , Ai < a < Aq 
= o V(u A ,A)er A . 



5. 



lim 

A— A* 



l ll|Y 



(2) There is a neighborhood U C X of u = 0, such that for each A with Ai < 
A < Aq + 6 and A ^ Ag, U contains only two nontrivial singular points n A 
and u\ of ( A. 16 ). 

(3) For each Ai < A < Aq + <5, U can be decomposed into two open sets U = 



U£ + U£ with n f/ 2 A = 0, such that 
(a) ifXi < A < A§, 

e t/ x x , 



€ J7. 



2 • 



t A e bu£ n <9[/ A , 



wi£/i u = and n 2 6emj attractors which attract t/ A and U2 respec- 



tively, and 
(b) ifX s <X<X s + 6, 

4 e (7x A , 



«2 e E# , 



e <9c/i A n dt/ 2 \ 



wi£/i n A and u 2 oeing attractors which attract /7 A and 1/^ respectively. 
(4) Near (u, A) = (0, Aq),u a and u 2 can 6e expressed as 



(A.20) 



i A = ai(A,e)e + o(|ai|), 
a 2 (A,e)e + o(|a 2 |), 



«2 



lim «i(A, e) = 0, 
A^Ag 

a 2 (Ao,e) ^ 0, 



where e is as in |j4.i9). 



Theorem A. 9. Assume the conditions (A. 4) and (A. 5) with m = 1, and (A.19) 

with b\ > 0. Then, there is an e > s«c/i £/ia£ when and T| satisfy (A. 17), 
the transition of (A. 16) is either jump or mixed. If it is jump transition, then 
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Assertions (1) and (3) of Theorem A. 4 are valid for (A. 16). If it is mixed, then 
the following assertions hold true: 



(1) (A. 16) has a saddle-node bifurcation at some point {u\, \\) € Xx (A§, +oo), 

(i = l,2), 



and there are exactly two branches 



Tt = {(«?, A)| \^-6<\< Ai} 
separated from (ui,Ai), which satisfy 

\\v£\\x = o V(n A ,A)eL*, 
o V(n A ,A)eL A . 



AS — e < A < Ai, 



lim || u 

A^Ag 



A ll 

i IIa 



(2) There is a neighborhood U C X ofu — Q, such that for each A with Xn__6< 
A < Ai, U contains only two nontrivial singular points n A and u\ of (A. 16). 

(3) For every Ag — 5 < A < \±,U can be decomposed into three open sets 
U = U Q + U x + U 2 with Ui f~l Uj = (i ^ j) such that 

(a) if Aq — (5 < A < Aq, £/ien 

« = oe c/ \ u A e <9[/ A n <9£/ A (i = 1,2), 

urai/i u = feeing an attractor which attracts Uq and U*(i — 1,2) iwo 
saddle points with the Morse index one, and 

(b) if Aq < A < Ai, iften 

n A e f7 x A , u A e <9t/ 2 A n flc/^, o e 8U} n ac#, 

u A being an attractor which attracts £/ A and u A and n = being 
saddle points with the Morse index one. 



(4) Near (0, A°),u A and u A can oe expressed by (A. 20) 



Appendix B. Ginzburg-Landau Models 

In this subsection, we introduce the time-dependent Ginzburg-Landau model for 
equilibrium phase transitions. 

We start with thermodynamic potentials and the Ginzburg-Landau free energy. 
As we know, four thermodynamic potentials- internal energy, the enthalpy, the 
Hclmholtz free energy and the Gibbs free energy-are useful in the chemical ther- 
modynamics of reactions and non-cyclic processes. 

Consider a thermal system, its order parameter u changes in f2 C K n (1 < n < 3). 
In this situation, the free energy of this system is of the form 



(B.l) 



H(u,X) =H 



1 m 



/iilVt^l 2 +g(u,Vu,X) 



dx 



where N > 3 is an integer, u — (ui, • • • , u m ), fa = fa(X) > 0, and g(u, Vn, A) is a 
C"*(r > 2) function of [u, Vu) with the Taylor expansion 



N 



(B.2) 



g(u, Vu, A) = ^2 ®ijkUiDjU k + aju 1 + o(\u\ N ) - fX, 

|/|=i 



where / = (h,--- ,i m ),ik > are integer, |/| = X)l'Li*fci the coefficients a^fe 
and a j continuously depend on A, which are determined by the concrete physical 



problem, u 1 — w^ 1 



and fX the generalized work. 
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A thermal system is controlled by some parameter A. When A is far from the 
critical point Ao the system lies on a stable equilibrium state £1, and when A 
reaches or exceeds Ao the state Si becomes unstable, and meanwhile the system 
will undergo a transition from Si to another stable state S2. The basic principle 
is that there often exists fluctuations in the system leading to a deviation from the 
equilibrium states, and the phase transition process is a dynamical behavior, which 
should be described by a time-dependent equation. 

To derive a general time-dependent model, first we recall that the classical le 
Chatelier principle amounts to saying that for a stable equilibrium state of a system 
S, when the system deviates from S by a small perturbation or fluctuation, there 
will be a resuming force to restore this system to return to the stable state S. 
Second, we know that a stable equilibrium state of a thermal system must be the 
minimal value point of the thermodynamic potential. 

By the mathematical characterization of gradient systems and the le Chatelier 
principle, for a system with thermodynamic potential 7i{u, A), the governing equa- 
tions are essentially determined by the functional Tl(u, A). When the order param- 
eters (tti, • • • , u m ) are nonconserved variables, i.e., the integers 

Ui(x,t)dx = a,i(t) ^ constant. 

n 

then the time-dependent equations are given by 
du 
~di 

(B.3) du 



= -0ij-H(u, A) + ^(u, Vu, A), 
at ouj 



n ,af2=0 (oru|ao = 0), 
u(x, 0) = ip(x), 

for any 1 < i < m, where 5/5ui are the variational derivative, Pi > and satisfy 

(B.4) / V H(u,X)dx = 0. 

Jn , oui 



The condition (B.4 1 is required by the Le Chatelier principle. In the concrete 



problem, the terms <&i can be determined by physical laws and (B.4). 

When the order parameters are the number density and the system has no ma- 
terial exchange with the external, then uj (1 < j < m) are conserved, i.e., 

(B.5) / Uj(x, t)dx = constant. 

Jn 

This conservation law requires a continuous equation 

nit 

(B.6) ^ = _ V .J J ( U ,A) > 

where Jj(u, A) is the flux of component Uj. In addition, Jj satisfy 

(B.7) J^-kjVfa-^fn), 

where /ii is the chemical potential of component ui, 

(B.8) fij -^2m= —H(u, A) - <j)j(u, Vm, A), 
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0. 



and 4>j(u, A) is a function depending on the other components Ui (i ^ j). When 
m = 1, i.e., the system consists of two components A and B, this term 
Thus, from (B.6l-(B.8l we obtain the dynamical equations as follows 

8 



du. 



(B.9) 



■ = PA 

an = 0, 



dAu, 



at 

du 

on on 
u{x, 0) = <p(x), 



TL{u, A) — ^ (m, Vm, A) 



I an 



0. 



for 1 < j < m, where ft > are constants, ^ satisfy 



(B.10) 



Jn ., 



5w 



-H(u,X)dx = 0. 



If the order parameters (ui,-- - ,Uk) are coupled to the conserved variables 
(uk+i, • • • , u m ), then the dynamical equations are 

= -ft j-H(u, A) + (u, Vu, A) , 



(B.ll) 



9uj 

at 

dn 
duj_ 
dn 



PA 



-7i(u, A) — <pj(u, Vu, A) 



I an 



\en 



0, 



(or Ui|gn 



9n 



<«■> 



0), 
= 0, 



u(x, 0) = if(x). 

for 1 < i < k and fc + 1 < j < m. 

The model (B.ll I gives a general form of the governing equations to thermody- 
namic phase transitions. Hence, the dynamics of equilibrium phase transition in 
statistic physics is based on the new Ginzburg-Landau formulation (B.ll I. 

Physically, the initial value condition it(0) = (f in (B.ll) stands for the fluctu- 
ation of system or perturbation from the external. Hence, <p is generally small. 
However, we can not exclude the possibility of a bigger noise ip. 
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